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Q\ '. Abstract 



\^ • We prove a weak law of large numbers for a tagged particle in a 

^P . totally asymmetric exclusion process on the one-dimensional lattice. 

The particles are allowed to take long jumps but not pass each other. 
The object of the paper is to illustrate a special technique for proving 
such theorems. The method uses a coupling that mimics the Hopf- 
Lax formula from the theory of viscosity solutions of Hamilton-Jacobi 
equations. 
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1 Introduction 

The purpose of this paper is to explain a technique for deriving scahng hmits 
that apphes to certain totally asymmetric particle processes. The example 
we use to illustrate the method is a one- dimensional exclusion process whose 
particles are allowed to take arbitrarily long jumps, as long as particles do 
not pass each other. Thus the ordering of particles is preserved, and the 
process models traffic on a single-lane highway. To the author's knowledge, 
the hydrodynamics of this particular process have not been addressed before. 

The key point of the technique is this: There exists a coupling that ex- 
presses a process with general initial conditions in terms of a family of pro- 
cesses that have a simple structure. The evolution from a general initial 
configuration is constructed from the simple cases by a variational formula, 
hence the term "variational coupling" . The benefit of the scheme is that the 
limit behavior of the simple processes is tractable, and through the coupling 
one obtains results for the general process. 

This paper is intended as an introduction to the technique. Fairly com- 
plete proofs are provided, and we have not strived to obtain the best possible 
result for the process under study. Some extensions of the theorem of the pa- 
per and shortcomings of the method are addressed in Section |^. A review of 
limit theorems for tagged particles and references to literature can be found 
ini. 

Notational remarks. For a real number x, [x\ denotes the largest integer 
less than or equal to x. Z_ is the set of nonpositive integers {. . . , —3, —2, —1, 0}. 

2 Description and construction of the process 

The process consists of particles that move on the one-dimensional lattice 
Z. At most one particle is permitted on each site. (This is the meaning 
of the name "exclusion process".) We label the particles by integers i G Z. 
The state of the process is a sequence of integers a = {cr{i) : i G Z), where 
a{i) is the location of particle i. We shall also use "cr(i)" as the name of the 
particle, because typically we need to discuss simultaneously several different 
processes denoted by different Greek letters. The dynamics will be such that 
the ordering of particles is preserved. So we can assume that 

a{i) + 1 < a(i + 1) for all i G Z. (1) 



This condition contains both the exclusion rule and the no-passing rule. 

It is also convenient to allow the values ±00 for cr{i). Assumption ([^) 
remains valid. If cr{i) < cxo = a{i + 1), then a{i) is the rightmost particle on 
the lattice Z because a{j) = 00 for all j > i. Similarly there is a leftmost 
particle if some cr(j) = —00. 

To summarize, the state space of the process is the set of ZU{ ±00}- valued 
sequences a that satisfy (|lD. 

For the dynamics, assume given a sequence Pk, k > 1, oi nonnegative 
real numbers. These are the rates. For each k and each i, particle a{i) 
independently attempts to jump k steps to the right at exponential rate Pk- 
If there are k empty lattice spaces in front of cr(z), this jump can be executed. 
If there are not, in other words, if a{i + 1) < <7{i) + k, then a{i) jumps as 
much as it can, to the location a{i + 1) — 1 right behind the next particle. 
The jumps occur independently for all particles a{i) and all jump sizes k. 

For nonnegative integers m, define 

00 

Bm = T. ^"A ■ (2) 

fc=l 

We assume 

B2 < 00. (3) 

Then also -Bq < 00, which guarantees that there are finitely many jump 
attempts in finite time intervals almost surely. B2 < 00 implies that the 
distance traveled by any particle during a finite time interval has a finite 
second moment. 

To rigorously construct the process, let {Vi : i e Z} be an i.i.d. collection 
of homogeneous rate Bq Poisson point processes on [0, cxo). In other words, 
a time interval [t, t + h) has a Poisson(i?o^)-distributed random number of 
points (or epochs) from the process Pj. For each process Vi, label the points 
independently by integers k > 1, with probability (3k/ Bq for label k. Let T)^ 
denote the point process of fc-epochs (that is, the points with label k) from 
process Vi. It follows that {V^ : i e Z, /c > 1} is a collection of mutually 
independent homogeneous Poisson point processes on [0, 00), and the rate of 
V^ is Pk- If /5fc = the point process Pf is empty for all i. 

The idea of the construction is that particle a{i) reads its jump commands 
from the Poisson process Vi, and attempts to execute a fc-step jump at each 
epoch of Pf . Let a{i, t) denote the location of particle a{i) at time t > 0, and 



a{t) = {(j{i, t) : 2 G Z) the entire state at time t. Assume that a deterministic 
initial configuration cr(0) = {(T{i, 0) : i G Z) has been specified. Informally, 
the rule of evolution is this: 

Jump rule. Suppose r is an epoch of Vf. Assume that the state a(t) has 
been defined for t < t. At time r set 

cr(i, r) = min {o"(i, r— ) + k , a{i + 1, r— ) — 1} . (4) 

After r, a{i) stays constant until the next epoch in Pj. The same scheme 
happens simultaneously and independently for all particles cr{i). 

Since there are infinitely many particles, infinitely many jumps are at- 
tempted almost surely in any time interval (0, e). In particular, there cannot 
be a first jump. Consequently some argument is needed to justify the induc- 
tion in the jump rule, so that formula (|D can be used to define the successive 
locations of particle cr(i). 

The realizations {T>^} for which we can construct the dynamics satisfy 
these assumptions: 

(a) Each Vi has only finitely many epochs in each bounded 
time interval. 

(b) There are no simultaneous jump attempts. 

(c) There are arbitrarily large io and to such that V^i^ and Vi^ 
have no epochs in [0,to]- (5) 

It is a standard fact of Poisson processes that these assumptions are satisfied 
by almost every realization {"Df}. 

Assumptions @ allow us to prove that the jump rule defines the evolution 
unambiguously for all particles a{i) and all times < t < oo: Given i and t, 
part (c) of d^) gives io and to so that —i^ < i < iq and tg > t, and so that 
P_Jo ^'^d PjQ have no epochs in [0,to]. Consequently cr(— 2o,t) = cr(— io,0) 
and cr(2o,t) = cr(2o,0) for all t G [0,to], and the evolution of particles a{i) 
for — io < j < io is isolated from the rest of the process up to time to. By 
parts (a) and (b) of (||), the locations cr(j, t) for —i^ < j < io and t G [0,to] 
can be computed by applying the jump rule to the finitely many potential 



jump times in U-ig<:j<:igVj fl [0, to], in their temporal order. In particular, the 
motion of particle a{i) is defined up to time t. 

This way the evolution a{-) = {<j(t) : t > 0} is defined as a function of 
an arbitrary initial configuration cr(0), and of almost every realization of the 
jump times {V^}. If (f2,jF, P) is a probability space on which are defined 
the point processes {V^} and a random initial configuration a"(0), then o"(-) 
is a stochastic process defined on this same probability space. Note: The 
Poisson point processes have to be independent of the initial configuration 
a(0). 

Constructing particle systems from Poisson processes of potential jump 
times is a standard procedure, and can be found in general references on 
particle systems such as [|I], |, |^. A useful property of this construction is 
that it preserves orderings between two processes: 

Lemma 1 Suppose the probability space {fl, JF, P) contains two initial con- 
figurations cr'(- ,0) and a"{- ,0) that satisfy cr'(z,0) < cr"(i,0) for all i E Z 
with probability 1. // the processes a'{- ,t) and a"{- ,t) obey the same point 
processes {T>^}, we have cr'{i,t) < a"{i,t) for all i E Z and t > 0, with 
probability 1. 

Proof Choose io and to to satisfy part (c) of assumption (^). We shall 
prove that the ordering a'{i,t) < a"{i,t) holds for —io < i < io and t < to. 
It holds for i = ±io because particles a'{±io) and a"{±io) do not attempt to 
jump during the time interval [0,to]. By parts (a) and (b) of (H), we may do 
induction over the jump times. Suppose r G Pf is the first epoch at which 
the ordering is violated, for —io < i < io and r < to. Then we must have 
a'{i,T) > a"{i,T) but a'{i,T—) < a"{i,r—). Since cr"(i) was not allowed to 
jump as far as a'{i) even though they both attempted a fc-jump, it must be 
that a'{i + 1,t— ) > a"{i + l,r— ). Since r is the first time the ordering is 
violated for indices — io < "^ < ^o, it must be that i + 1 > i^. Since io > ^5 
we conclude that i = io- But Pj^ was chosen to have no epochs up to time 
to. We have reached a contradiction, so there can be no first violation of the 
ordering for — io < i < io and t < to- Since io and to can be taken arbitrarily 
large, the proof is complete. I 



3 The heuristic behind the method 

We are interested in a macroscopic view of the process. By this we mean 
its large scale behavior, distinct from the local interactions of the individual 
particles a{i) that represent the microscopic dynamics. The ratio between 
macroscopic and microscopic scales is given by a parameter n, so that macro- 
scopic space and time units represent n microscopic space and time units. We 
hope that a deterministic macroscopic description of the dynamics emerges 
in the limit n ^ oo, where the macroscopic and microscopic scales become 
infinitely separated. 

What form should the deterministic macroscopic description take? Sup- 
pose we have a function u{x,t) defined for macroscopic space-time variables 
(x, t) such that the approximate equality 

u{x,t) ^ n^^a([nx],nt) (6) 

holds with high probability for large n. Since the rate of advance of particle 
a{[nx],nt) depends only on the distance to the next particle on the right, we 
might expect this same structure also at the macroscopic level, so that 

ut = fM (7) 

for some function /. Equation (^ is a first-order partial differential equation 
of the Hamilton- Jacobi type. It is usually written as Ut — f{ux) = 0, and 
then the function —/is known as the Hamiltonian. 

The exclusion rule and (|^) force m^: > 1, so / is a function defined on 
[1, oo). If the distance to the next particle ahead is one (no empty sites) a 
particle cannot jump, and hence /(I) = 0. The rate of advancement cannot 
decrease as the distance to the next particle increases, so / ought to be 
nondecreasing. No particle can advance faster than at rate Bi, so f < Bi. 
Since / is nondecreasing and bounded above, it appears reasonable to assume 
that / is concave. 

Thus, assuming our heuristic reasoning is reliable, if there is to be a law 
of large numbers 

n~^a{[nx],nt) -^ u{x,t) , (8) 

then u should satisfy equation (0) with a convex Hamiltonian — /. There is 
a well-known formula for solving (^ when the Hamiltonian is convex: Let 

g{x) = sup{fx + /(f)} (9) 

V>1 



be the convex conjugate of — /. Suppose the initial data for (|^ is given by 
a function uq on R. Then the Hopf-Lax formula 

uix, t) = inf i^uoiy) + tg (^) } (10) 

defines the unique viscosity solution of (0) with initial data u{x, 0) = uo{x) 
(Theorem 3, Section 10.3 in 0). Note: Since / > 0, (H) shows that g{x) = cx) 
for X > 0. That is why we may restrict the infimum in (|To|) to y > x. 

However, our task is not to solve the differential equation (^ with a given 
/. In some sense we already have the solutions, at least in their microscopic 
form a{[nx],nt) because we already constructed the particle process. Instead, 
our problem is to prove that the process a admits a macrosopic description of 
the type (|^), and to find / and g, if possible. The Hopf-Lax formula suggests 
a line of attack on our problem. 

Suppose we are allowed to feed a single input uq = i^q to the system 
described by (p!OD, and measure the output (p{x,t) = u{x,t). Can we obtain 
g7 Take 

^^^y^ = [y if'<o: (11) 

With Mo = v^o; differentiating inside the braces in ([T0| ) shows that ^{x,t) = 
u{x,t) = tg{x/t). (The duality (|) implies that the slope of g is everywhere 
at least 1, since this is the lower bound for v.) Once we have g, we can 
compute the solutions to (^ by the Hopf-Lax formula (pJOf ). In other words, 
the single evolution (p{x,t) started with ([TT]) contains the information for 
constructing all evolutions from their initial data. 

To understand this point better, let mq be an arbitrary initial function 
for the evolution, and consider this family of special initial data, indexed by 
z e R: 

if'oiy) = Uo{z) + ipoiy) , yen. (12) 



These are simply translates of the initial data (|TlD. Find the evolution from 
(HOD: 



ip'ix,t) = Mi^^',{y)+ig[^j^) 



Uo{z) + mi \ipo{y)+tg 

y>x L \ t 

Uoiz) +tg{x/t). 



x-y 



We substitute the evolution y9^(x, t) into (|T0|) to obtain 

u{x, t) = inf ip'^{x — y,t) . (13) 

y>x 

Equations ([T^) and ([T3| ) contain the message: A general evolution u{x,t) 
starting from initial data uq can be computed from the family of solutions 
{ip^ : y G R}, whose initial data are translates ip^ = Uo{y)+ipo of the function 

Let us move the discussion to the microscopic particle level. We can 
realize the initial data (|ll]) microscopically, in the sense of the limit (H), by 
a natural particle configuration: 

.. ^. f oo if i > 0, ,, .. 

^i^,0) = [, if,<o: (14) 

This says that the lattice is empty to the right of the origin, while to the left 
all sites are occupied, the origin itself included. Since (|1^ is a microscopic 
version of ([TT|), we would expect that the process a{t) with initial configura- 
tion (|T^ is a microscopic version of the solution ip{x, t) = tg{x/t). Translates 
of these solutions suffice to express the most general solution, as observed in 



The key question becomes: Is there a microscopic principle that corre- 
sponds to ( |T3|) and expresses a process with a general initial configuration 



in terms of a family of processes of the type that start from (0)? The 
construction we seek is the coupling we now turn to. 

4 The coupling 

Assume given an arbitrary deterministic initial particle configuration a{0) = 
(cr(j, 0) : j G Z), and let a{j, t) denote the process constructed in terms of the 
Poisson point processes {T>^} as explained in Section ^ For each finite initial 
location 0"(j, 0) we construct an auxiliary exclusion process C,^ = {(^{i,t) : 
i G Z_,t > 0). Initially 

(j(i^ 0) = a{j, 0) + i for i G Z_. (15) 

The particles C"'(0 ^^ ^ > are not needed. According to our convention we 
may think that they reside permanently at oo. Thus (^ is a process that is 
a spatial translate of the process started from (p!^. 



For each process (^ , the dynamical description is the same as for a: Jumps 
of size k are attempted independently at rate /3fc. Each jump is carried as 
far as possible without violating the exclusion and no-passing rules. Each 
process (^ uses the same collection of Poisson processes {^f } as does a, but 
with a translation of the index: 

Jump rule for process (^ . At each epoch r oi'D'^j^-, set 

C^{i, t) = min [C'{i, r-) + k, CK^ + 1> r-) - l} . (16) 

Between epochs of Pj+j, CK'') stays constant. 

Subject to this rule, the processes {C,^} are constructed exactly as a. The 
processes a and {Q} are all defined on the same probability space (f2, JF, P) 
of the Poisson jump times. If this probability space also supports a random 
initial configuration cr(0), then the processes {C,^ are functions of both o"(0) 



and {^j }, through the translations (|T^) and the jump rule 



This arrangement where several processes are defined on a common prob- 
ability space to facilitate their direct comparison is known as a coupling. It 
should be evident that the processes are invisible to each other, i.e. particles 
of one process do not interfere with the jumps of the other processes. 

The key property of this coupling is the microscopic counterpart of ([13|) : 

Lemma 2 The equality 

aii,t)= mlC{t-j,t) (17) 

holds for all i E Z and t > 0, almost surely. 

Proof. The exclusion rule (|I|) and (|T3p imply that (|T^) holds at time 0. 
The jump rules ensure that for each j > i, particles C'^l^^j) and a{i) attempt 
fc-jumps at common time points, namely at the epochs of Pf . Thus for a 
fixed i the validity of (|^) can change only at epochs of Vi. 



Fix ii and ti. To prove that ([17D holds for i = ii up to time ti, use part 
(c) of assumption (||) to find io such that —io < ii < io, and so that P-jp 
and Vig have no epochs in the time interval [0,ti]. By parts (a) and (b) of 
(^, we can do induction on the finitely many epochs in U_i„<j<j(,'Dj fl [0,ti]. 

9 



So suppose r is an epoch of V^ for some A; > 1 and some i between —io and 



io, and assume that ([T7|) holds for all i between — io and io and all t < t. 

First we show that < holds in ([T^ right after the jump epoch r. Let j > i 
be arbitrary. By the induction assumption, (^{i — j + 1, r— ) > a{i + 1, r— ), 
and consequently (^{i — j) can jump at least as far as a{i) can. Since before 
the jump we have C"'(i ~ j,T—) > a{i,T—), and both particles attempt a 
/c-jump, the inequality (^{i — j,r) > cr(i,r) holds also at time r after the 
jump has been completed. The same argument works for all j > i. 

Next we show that right after r there is some j > i that satisfies (^{i — 
j, t) = a{i, t). If cr(i) was allowed to jump the full k steps, any j such that 
C^^ii — j, r— ) = a{i, r— ) will do, and by induction there is some such j. 

Suppose on the contrary that a{i) was at least partially blocked by o"(i+l) 
at time r— , so after the jump we have a{i, r) = a{i + 1, r) — 1. By induction, 
we can pick j > i + 1 such that C,\i—j + l^ t—) = cr(i + l, r— ). By the earlier 
part of the proof, (^{i — j, r) > cr(i,r). On the other hand, the exclusion 
rule for the process (^ stipulates that (^{i — j, r) < (^ {i — j + 1, t) — 1 = 
a{i + l,r) — 1 = a{i,T). (Notice that by part (b) of assumption (^ there is 
no epoch r in Pj+i.) This shows that after the jump we necessarily have the 
equality (^ (i - j^r) = a{i,T). I 

By subtracting off the centering at o"(j, 0) from C,^ we obtain a family of 
processes {^^} that are all equal in distribution. Define 

i^{i, t) = CKh t) - crU, 0) for all t,j, and t. (18) 

Then we can turn the coupling equality ([T7| ) into a statement that is a mi- 
croscopic version of the Hopf-Lax formula (PUJ): 

a(z,t)= inf.{a(j,0)+e(z-j,t)}. (19) 

The processes {^■'} are identically distributed on account of the requirement 
that (o"(j, 0)) and the Poisson processes {Vj} are independent. The advan- 
tage of ( [T9| ) over ( p!7D is that the effects of (cr(j, 0)) and {Vj} have been 
separated in two distinct terms inside the braces. 

5 The tagged particle limit 

As a corollary of the coupling we obtain the limit discussed heuristically in 
Section ^. The precise assumptions are as follows: There is a sequence of 

10 



processes an constructed as in Section ^ and an increasing function uq on 
R. At time 0, the convergence 

\imn-^ani[ny],0)=Uoiy) (20) 

holds in probabihty for each y E H. More precisely, if (f2„,jF„,P„) is the 
probability space on which the process (j„ is defined, then 

limPn{\an{[ny],0)-nuo{y)\>ne) = (21) 

for each y G R and e > 0. 

Theorem 1 Under assumption (pO]), the convergence 



Mm. n a„i\nx],nt) = u(x,t) (22) 

/io/ds m probability for each a; e R and t > 0. The deterministic limit u{x, t) 
can be macro scopically defined by 

u{x, t) = mi^ ^u^{y) + tg (^^) } , (23) 

where g is a certain convex function on (— oo, 0] determined by the rules of 
the exclusion process. 



The meaning of the limit in (|22| ) is the same as above for (pO[), namely 
that 

lim Pn (|(T„([nx],nt) — nu{x,t)\ > ne) = (24) 

for each e > 0. By the discussion about Hamilton- Jacobi equations in Section 
1^, u{x,t) can be equivalently characterized as the unique viscosity solution 
of Ut = f{ux) with initial data ^(x, 0) = uo{x), where / is the negative of 
the convex conjugate of g. The theorem says that in this exclusion system, 
a tagged particle has a well-defined macroscopic velocity f{ux) as long as 
initially the particle distribution is sufficiently regular as specified by as- 
sumption (pop. 



11 



Proof of Theorem 1 

To prove Theorem |I|, we forge a rigorous connection between the matching 
microscopic and macroscopic descriptions (|l^) and (^31). For each n, we have 
the processes C,^ and ^^, defined as in Section ^ relative to the process an- 
All the processes ^ are identical in distribution, so as long as we discuss 
distributional properties, we may leave out the sub- and superscript n and 
j. The function g that appears in the theorem is defined by the limit of the 
next lemma. 

Lemma 3 There exists a finite, continuous, convex function g on (— cxd, 0] 
such that 

lim n'^^{[nx],nt) = tg{x/t) (25) 

in probability, for each x < and t > 0. 

Proof of Lemma ||. We start by deriving a useful subadditivity. Assume 
the process ^{i,t) has been constructed on a probability space (f2, JF, P) as 
explained in Section |^, so that particle ^{i) gets its jump commands from the 
labeled Poisson process T>i. To perform comparisons, we explicitly include 
the particles that are permanently at infinity: C,{i,t) = cxd for alH > and 
t >0. 

Fix an integer /i G Z_ and a time s > 0. Construct two new processes 
a' and a" on the same probability space {Q, JF, P) as follows: The initial 
configurations are 

a'{i, 0) = ^{h + i, s) for alH G Z (26) 

and 

^>^'^ = {ah,s)+^ if.<0: (27) 

We stipulate that particles cr'{i) and cr"{i) read their jump commands from 
{Vh+i — s) n [0, oo). In other words, a' and a" read the Poisson processes 
{'Di} after first translating the index by h and time by s. By Lemma |I|, 

a'{i, t) < a"{i, t) for alH G Z and t > 0. (28) 

Since a' continues ^ with translated index and time, a'{i, t) = ^{h + i,s + t). 
For a" we can write cr"{i, t) = ^{h, s) +C,{i, t) where ^ is equal in distribution 

12 



to C, but independent of it. Thus ( pH] ) turns into 

ah + t,s + t)<ah,s) + i{t,t), (29) 

where the random variables on the right are independent. 

Suppose first that x in ( ]25|) is a nonpositive integer, x = i. Let F„ be 
the distribution function of ^{ni,nt). Then (|29|) gives -F^+m > -F'n * -^m- By 
the Kesten-Hammersley lemma from subadditive ergodic theory there exists 
a function 7(^,t) such that 

lim n~"^^(m, rat) = 7(z,t) (30) 

in probability, for all i G Z_ and t > 0. A proof of this lemma can be found 
on p. 20 in |jll|. The assumption of finite second moments required by this 



lemma is satisfied by assumption (|). 

The rest of the proof consists of extending the limit to all x < 0, by using 
regularity and monotonicity properties of the process and 7(«, t) . It is obvious 
that 7(z, t) is nondecreasing in both i and t. For ? = we have 7(0, t) = Bit, 
because Bi is the speed of an unobstructed particle. Furthermore 

i<-fii,t) <i + Bit (31) 

for all i G Z and t > 0, because C,{i, 0) = i and Bi is the maximal average 
speed. 

Next we argue the t-continuity of 7(i, t). During the time interval (nt, nt+ 
ne] particle ^{ni) attempts a Poisson(i?on£:) number of jumps. The attempted 
jump sizes are i.i.d. random variables {Xi} with common distribution P(Xj = 
k) = 13k/ Bq. Due to the exclusion and no-passing rules, the actual jumps are 
at most Xi, so we can write 

N 

^{ni,nt + n€)<^{ni,nt) + ^X^ (32) 

where A^ ~ Poisson(5o?^£). It follows that 

'y{i,t + e)<-f{i,t) + Bie. (33) 

In particular, 'y{i,t) is Lipschitz continuous in t. 
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A homogeneity holds: If m is a positive integer, then 

'y(mi, mt) = hm n~^f(nmi, nmt) 

= m ■ liin (n'm)~^f(n'mi, nmt) 
= m^ii^t) . 

From this it follows that we can define unambiguously for rational r < 

7(r, t) = — ■y(mr,mt) , (34) 

m 

where m is any positive integer such that mr G Z_. 

We need to prove that the extension (|3^) makes sense for the process, 
namely that 

lira n~^f(\nr],nt) = 'J (r,t) (35) 

n^oo 

for rational r < and all t > 0. If n runs along the subsequence n = km, 
k = 1,2,3, ... , where m is the integer in (0), then (|35|) is valid. For other 



values of n we interpolate using the monotonicity of the process. No-passing 
and total asymmetry of the dynamics imply that 

for all i < j and s < t. Given n, let k = k{n) be the unique integer that 
satisfies km, <n<{k + l)m. Let £ > 0. Then for large enough n, 

^{{k + l)mr,{k + l)m{t-£)) < i{\nr],nt) 

< C, {kmr, km{t + e)) , 

from which, upon dividing by n and passing to the limit, 

— 7 (?7ir, m(t — £:)) < liminf n~"'^^([?2r],nt) 

< lim sup n~^^([nr],nt) 

n— »oo 

< — 7 (?7ir, r7i(t + e)) . 
m 

These inequalities cannot be taken literally since the limits have to be inter- 
preted in the sense of convergence in probability. But the point should be 



clear: By the continuity (|33| ) and the definition (|34D, the limit ( pSf ) is valid. 
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Next one checks that the homogeneity extends to rationals: 

7(rx, rt) = r7(x, t) (36) 



for rational x < 0, rational r > 0, and all t > 0. From p9| ) follows subaddi- 

tivity 

l{x + y,s + t) <-f{x,s)+-f{y,t), (37) 



first for integers x, y and then by ( plf ) for rational x, y. (^) and (|37|) together 
imply convexity, and also continuity in the x-variable, for rational x. 
The final extension of 7 is 

7(x, t) = sup{7(r, t) : r G Q, r < x} . (38) 

By the continuity in the x-variable, for rational x this definition gives the 
old value 7(x, t), so it is a sensible extension. Again, one checks that homo- 



geneity ([36| ) and subadditivity (|37D hold, this time for all r, x, y. Convexity 
follows, and convexity implies continuity in the open quadrant {a; < 0, t > 0} 
(Theorem 10.1 in 0). Continuity up to the boundary requires separate ar- 
guments. And again, one checks that the limit ( ^5]) holds also for irrational 
r. This follows as it did for rational r, by the monotonicity of the process 
and the continuity of the limit. 

Finally, define g{x) = 7(x, 1). By homogeneity 7(0;, t) = tg{x/t). I 

We return to the proof of Theorem [^. Fix (x, t). From the coupling (|T9|) 
we get 

n^^an{[nx],nt) = inf ln^^an{i,0) + n^^^l^{[nx] — i,nt)\ . (39) 

i:i>[nx] ^ J 

By the hypothesis of Theorem |l| and by Lemma ^, for i = [ny] the random 
variable inside the braces converges to the quantity inside the braces in (|23|) . 
Consequently we get 

lim P„ ('n~V„([nx],nt) < u{x,t) + e) = 1 (40) 

for any e > 0. 

For the converse we need some further estimation. The first step is to 
restrict the variable i in (|39D to a range of order n. This is achieved by the 
next lemma. Define 

Wniz) = min |(T„(z, 0) + ^^ ([nx] - z, rat)} . (41) 

i:lnx]<i<[nz] ^ -' 
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Lemma 4 For large enough z > x, there exists a constant C = C{z) > 
such that 

Pn {an{[nxlnt) ^ Wn{z)) < e"^" (42) 

for all large enough n. 

Proof. Suppose particle ^jj"^' {[nx] — [nz]) has not moved by time nt, in 
other words that ^l^"^! {[nx] — [nz].,nt) = [nx] — [nz]. Let i > [nz]. By the 
exclusion rule for a^ and because 

Cn (M - h nt) > Cn i[nx] -i,0) = [nx] - i , 

we have 

a„(z,0) +Cn {[nx] ~i,nt) 
> (Jn {[nz], 0) + z — [nz] + [nx] — i 
= an {[nz], 0) + [nx] - [nz] 

The conclusion is that the terms i > [nz] cannot contribute to the infimum 

Thus to prove the lemma, it suffices to show that for large enough z, the 
probability that ^j^"^! {[nx] — [nz]) has moved by time nt is at most e"*^". This 
is clear again on account of the exclusion and no-passing rules: The random 
time when ^^['^l {[nx] — [nz]) ffist moves has the distribution of a sum of 
[nz] — [nx] + 1 i.i.d. exponential random variables with rate Bq, because after 
^(j) ffist moves, C,{j — 1) waits an Exp(i?o)-distributed independent time to 
make its ffist jump, and these times are added up over j = 0,-1, —2, . . ., 
[nx] — [nz]. It is sufficient to choose z so that {z — x)Bq^ > t + e for some 

£>0. I 



Fix z so that (^2]) holds, and let e > 0. Pick a partition 
x = yo<yi<y2< ■■■ <ym = z 
so that 



, X — yk+i \ I X ~ yk 



, . , , < 7 (43) 

t J \ t J - A ^ ^ 

for k = 0,1, . . . ,m — l. If we could restrict the minimum in (^) to the m + 1 



points i = [nyk], then passing to the limit inside or outside the minimum (now 
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over a fixed finite number of random variables) would make no difference. To 
control the error of this simplification, we need to control what happens inside 
the braces in (^) for [nyk] < i < [nyk+i]- This can be done with couplings: 

Lemma 5 For alln andi < jo < ji, the inequality C,-I^{i— jo, t) > ^^{i—ji,t) 
holds almost surely for all times t > 0. 

Proof. Fix n and jo < ji- Consider the processes 

a'{i, t) = e(« - Ji, t) and a"{i, t) = e(^ - Jo, t), ^ G Z. 

Both cr'(z) and <y"{i) read their jump commands from Pj. (Recall the general 
rule from Section^: particle ^-^{i) reads jump commands from T>i+j.) Initially 

aUO) = |^ . 'J';^:^' (44) 

i «-Jl if « < Jl, 

and 

a'%0) = l°° . '!*>* (46) 

[ t - Jo It t < Jo- 

By the assumption jo < ji, <^'{h^) ^ <^"{h^), and consequently by Lemma 
|l|, the inequality cr'{i,t) < a"{i,t) holds for all times t > 0. I 

Apply Lemma || with i = [nx], jo = i between [nyk] and [nyk+i], and 
ji = [nyk+i], to bound Wn{z) from below by 

Wn{z) = min min |cr„(z, 0) + ,^^ ([nx] — ij^t) [ 

0<k<m i:lnyk]<i<lnyk+i] ^ J 



min {an {[nyk], 0) + d"'^+il {[nx] - [ny^+i], nt)} 



> 

Q<k<rn 



Multiplied through by n ^, the last line converges in probability by assump- 
tion (20) and Lemma ^ Thus with probability tending to 1 as n ^ oo, 

1 ^ \ ^ • f / N , ^ fx-yk+i\\ e 



n ^Wn{z) > min \uo{yk)+tg 

0<k<in [_ 



t J ] 4 
x-yk\\ £ 
' 2 
X — y\\ e 

2 

e 

2 ■ 



= u{x,t) 
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This together with Lemma ^ gives 

hm Pn (n~^a{[nx],nt) > u{x,t) — £:j = 1 (46) 

for any e > 0. This completes the proof of Theorem |l|. 

6 Comments and extensions 

6.1 Almost sure convergence in Theorem [T| 

It may be possible to upgrade the convergence in Theorem |I] to almost sure 
convergence by developing large deviation estimates for the convergence in 



Lemma a. Such results can be found in IH1 and 10 . 



6.2 The role of invariant distributions 

Theorem 1 does not explicitly identify the function g or the macroscopic ve- 
locity /. This can be done when information about steady states is available. 
Consider the special case of the totally asymmetric simple exclusion process 
(TASEP) where /5i = 1 and j3k = ioi k > 1. Then it is known that if the 
interparticle distances are geometrically distributed with mean v G [1, cxo), 

P {a{i + 1) - a{i) = m) = (l - v'^y\'\ m = 1,2,3, ... , (47) 

an individual particle's jumps obey a Poisson process of rate 1 — v~^. From 
this information the function g can be calculated with the help of Theorem 
|1|, as we next show. 

Let each o"„ = a be a steady-state process with interparticle distances 
(47), and such that initially cr(0, 0) = 0. Then Theorem ^is valid for uo{y) = 



vy. Take x = and t = 1 in (pSJ) to get 

u{0, 1) = inUvy + g{-y)} = -g*{v) , (48) 

where g* denotes the convex conjugate oi g. Note: For the purposes of convex 
duality, functions should be lower semicontinuous and convex on all of R. 
(See 0.) The limit of Lemma ^ defines g on (— oo, 0], and we extend (? to R 
by defining g{x) = oo for x > 0. Then (ESj) is the same as 



g*{v) = sup{a;f — g{x)}. (49) 



Furthermore, for f < 1 we obtain g*{v) = oo directly from (^), by the 
observation that the slope of g is always at least 1. This is a consequence of 
the exclusion rule (|I|) and the hmit (^). 

Since a{0,n) is a Poisson((l — f~^)n)-variable, we can calculate explicitly 
that 

m(0, 1) = lim ra"V(0, ra) = 1 - v'^. 

Thus g*{v) = v~^ — 1 for v > 1. From this we obtain g by performing another 
convex conjugation: For x < 0, 

9{x) = g**{x) 

= sup{xf — g*{v)} 

V 

= supjxf — V^^ + 1} 
v>l 

= 1-2a/^. 



The tagged particle velocity is the function f{v) = 1—v^^ obtained above. 
The relationship of the mean interparticle distance v to the particle density 
p is p = v^^. As a function of density the velocity is v{p) = f{p^^) = 1 — p. 
From this we get the well-known macroscopic particle current of TASEP: 
pv{p) = p{l- p). 

TASEP was the process studied by Rost in 1981 in one of the seminal 
papers on hydrodynamic limits for asymmetric processes. With an approach 
completely different from ours he investigated the process C, and proved a 
limit theorem related to Lemma |^. The use of subadditivity is a common 
feature of both approaches. 

A thorough treatment of TASEP with the method of this paper appears 
in 0. In particular, when the invariant distributions are available for explicit 
calculations, in addition to the law of large numbers it is possible to obtain 
a (one-sided) explicit large deviation rate function for the tagged particle in 
a nonequilibrium process. 

6.3 More general jump rates 

We could admit the following additional feature to the process: Not only is 
a{i) blocked by the particle in front of him, but he is also held back by the 
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particle behind him, in this particular sense: For some constant L > 1, 

(r{i) <a{i-l) + L for all i e Z. (50) 

The proof of Theorem |1| works with suitable modifications. The initial con- 
figurations of the (■"'-processes have to be changed to (^{i) = o"(j, 0) + i for 
i < 0, and ("'(O = '^U^ 0) + -^^ ^^^ i > 0. This feature appears in |Q. 

Another generalization that this approach can handle is random rates. 
The reader may have noticed that very little depended on how the point 
processes Vf were actually produced, as long as they had some desirable 
properties such as those listed in assumption (^. In particular, the rates 
of the processes could vary so that the rate /3fc,i of Vf is a random variable 
whose realization is fixed. If this randomness is suitably systematic (i.i.d. 
for example) one would expect that some form of the subadditive ergodicity 
needed for Lemma ^ works again. Examples appear in [^] and [|l^ . 

However, so far it has not been possible to admit a more complicated 
dependence of the rate of jumping on the distance to the nearest particles. 
By a more complicated dependence we mean a general function f3{k, I, r) such 
that jumps of size k are attempted at rate (3{k,l,r) when there are / and r 
empty sites to the left and right. In the model of this paper this function is 
a step function: For a fixed k, j3{k, I, r) takes only two values depending on 
whether k < r or k > r. What appears to fail for more general rates is the 
proof of Lemma 0. 

Other shortcomings are that the method has been able to treat only 
totally asymmetric processes, and processes without passing. A more general 
asymmetric process would allow jumps both left and right, but with a higher 
rate for jumping to the right. 
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